arXiv:1509.02524v2 [astro-ph.CO] 31 Dec 2015 


Mon. Not. R. Astron. Soc. 000,[T](7]{2015) 


Printed 5 January 2016 (MN RTbX style file v2.2) 


The MeSsI (Merging Systems Identification) Algorithm & 
Catalogue. 


Martin de los Rios*, Mariano J. Dominguez R.*, Dante Paz, Manuel Merchan 

* Instituto de Astronomia Teorica y Experimental (CCT Cordoba, CON ICET, UNC), Laprida 854, X5000BGR, Cordoba, Argentina. 

^ Observatorio Astrondmico de Cordoba, Universidad Nacional de Cordoba, Laprida 854, X5000BGR, Cordoba, Argentina. 

^ Consejo Nacional de Investigaciones Cientificas y Tecnicas, Rivadavia 1917, C1033AAJ Buenos Aires, Argentina. 


Accepted XXX Received XXX; in original form XXX 


ABSTRACT 

Merging galaxy systems provide observational evidence of the existence of dark matter and 
constraints on its properties. Therefore, statisticaly uniform samples of merging systems 
would be a powerful tool for several studies. 

In this work we present a new methodology for the identihcation of merging systems and 
the results of its application to galaxy redshift surveys. 

We use as a starting point a mock catalogue of galaxy systems, identihed using FoF 
algorithms, which experienced a major merger as indicated by its merger tree. 

Applying machine learning techniques in this training sample, and using several features 
computed from the observable properties of galaxy members, it is possible to select galaxy 
groups with a high probability of having experienced a major merger. Next we apply a mixture 
of Gaussian technique on galaxy members in order to reconstruct the properties of the haloes 
involved in such merger. This methodology provides a highly reliable sample of merging 
systems with low contamination and precisely recovered properties. We apply our techniques 
to samples of galaxy systems obtained from SDSS-DR7, WINGS and HeCS. Our results 
recover previously known merging systems and provide several new candidates. We present 
their measured properties and discuss future analysis on current and forthcoming samples. 
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1 INTRODUCTION. 


Merging galaxy systems such as the Bullet Cluster dClowe et al. 
2006| >, Abell 520 ( [Mahdavi et al.|2007||Jee et al.|2012||Clowe et al. 


2012|JJee^^e£alJ20^, Baby BulletJBradac^^e£a]J200^, Pandora 
( [Merten et al.|2()ll7 7Musket Ball ( [Dawson et al.|201^ El Gordo 


([Menanteau et al.[[2012| [Ng et al.[[2014| [Molnar & Broadhurst 
2015[ l, Abell 1758 ( [Durret et al.|2005t and Abell 3716 jAndrade 


Santos et al.[[2015ll have provided observational evidence for the 


existence of dark matter. Most of them have been used to test 


the Cold Dark Matter (CDM) paradigm itself i 

Markevitch et al. 

([2004j, Hayashi & White 1 2006|l, Farrar & Rosen 

(2007 i,M 

ilosavl- 

jevic et al. ^2007[l, Springel & Farrar 1 2007 i,[Randall et al. 

i20^ 

Mastropietro & Burkert[ ( 2008k Lee & Komatsu[ j2010k 

Forero- 

Romero et al. 1 2010k Thompson & Nagamine 1 2012 1 , [Watson et al. 

( 2014k Thompson et al. (20 14^. Several statistical techniques have 


been proposed to measure dark matter properties using merging 


systems (such as the self interaction cross section. 

Massey et al. 

2011 

IHarvey et al.||2014| |Kahlhoefer et al.||2014| 

Harvey et al. 

2015 

1 . However, the lack of a complete and uniformly selected 


sample of merging systems prevents efforts to derive robust con- 


* E-mail:martin@oac.unc.edu.ar, mardom@oac.uncor.edu 


straints. In order to overcome this limitation different approaches 
have been proposed: the Merging Cluster Collaboration uses the 
radio emission due to induced shocks in the Intra Cluster Medium 
(ICM, [Feretti et al.|2012) to obtain high redshift merging system 
candidates. These systems have been studied using pan-chromatic 
observations and detailed merging kinematic Bayesian reconstruc¬ 
tions ( [Dawson|[2013[ l. X-ray imaging, spectra and Sunyaev Zel- 
dovich effect observations have been used to identify unrelaxed 
clusters of galaxies jMann & Ebelin g|2012[ l, cluster mergers ( [Har- [ 
[vey et al.|2014^ , substructures and any departures from hydrostatic 
equilibrium, mainly for the most massive galaxy clusters. 

Galaxy redshift surveys are very useful to trace the dynamical 
state of galaxy systems and to search for substructures (jDresslerj 
[& Shectman[[1988) l. To this end, some methods look for depar¬ 
tures in the global Gaussian redshift distribution of system mem¬ 
bers ( [Solanes et al.[[199^ [Hou et al.[[2014[ [Yu et al.[[2015p . Even 
though all these methods aim to identify the substructures and re¬ 
cover its properties, they suffer from false identifications and in¬ 
completeness, at least to some extent. 

In this work we develop a uniform identification algorithm of merg¬ 
ing systems based on galaxy redshift catalogues. These methods 
can be applied on low mass systems and should increase the num¬ 
ber of merging systems like the bullet group recently identified by 
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IGastaldello et al.| p014| l. This paper is organized as follows: In 
Sec|^ we apply machine learning techniques to a number of ob¬ 
servable features and present its calibration based on the result of 
simulations. We also introduce techniques to recover properties of 
merging dark matter haloes. In Sec.j^we apply our techniques to 
samples of galaxy systems identified from low-redshift galaxy sur¬ 
veys (SDSS-DR7 HeCS, WINGS). 

Finally, in Sec.|^we summarize the main results of this work 
and discuss uses of this new sample of merging systems. We adopt 
the standard cosmological model used in the Millenium Simula¬ 
tion jSpringel et al.|2005 1 when necessary {Hq = 73 kms^' Mpc^*, 
Sl„, = 0.25, and Oa = 0.75). 


members in order to have a better identification. We performed the 
DS test for « = 10 and n = and use global and individual 

DS values as features. 

(ii) Well known tests measuring the departures from a normal 
Gaussian distribution: The Anderson-Darlling test, the Cramer-von 
Mises test, Kolmogorov-Smimov test, Pearson chi-square test and 
Shapiro-Francia test (provided by the package nortest [Gross &] 
|Ligges|[2015^ and Shapiro-Wilk test (provided by the stats pack¬ 
age). 

(iii) Astrophysical properties of galaxies and clusters: SDSS 
magnitudes, g-r colour index and occupancy of clusters. 


2 METHODOLOGY. 

2.1 Mock Galaxy and Halo Catalogues 

From the point of view of the current theory of galaxy formation, 
the most direct route for defining galaxy systems is via its host dark 
matter halo ( |Mo et al.|2010^ . 

The Millennium simulation ( jSpringel et al.|2005| l, used in this 
work, provides a catalogue of dark matter (sub)haloes constructed 
using a traditional 3D Friends-of-Friends algorithm (FoF) that per¬ 
colates nearby particles. The GAVO Millennium data base, also 
provides merger trees for each halo ( jRoukema et al.|1993^ , there¬ 
fore given all the haloes belonging to a FoF group, it is possible to 
infer its merger tree. We compute for each FoF group a parent list 
of FoF progenitors (identified in previous snapshots) which have 
contributed with at least one subhalo to the current FoF group. We 
define a major merger of FoF groups as the merger between two 
groups where the total mass of the involved haloes represents at 
least a 20 per cent of the mass after the merger, it is worth noting 
that this condition imposes a minimal value of 0.25 for the mass ra¬ 
tio of the interacting systems. We construct mock catalogues of the 
SDSS-DR7 redshift survey based on the results of a semi-analytic 
model jGuo et JT]|2011[ ) and use it to calibrate our identification 
method for merging systems. This process is extensively described 
in previous works ([Lares et al.|201 IJ and (jDomlnguez Romero et al.j 

I20T2I 1. 

We define a recent major merger as a the major merger of two 
FoF groups, where its principal haloes are still present as different 
haloes in the final FoF group. With this selection criteria the mean 
of the look-back merger time is around 3 Gyr consistent with other 
work jPinkney et al.|[T996) . The FoF groups identified as recent 
major mergers and their member galaxies identified in the mock 
catalogs are used to train different machine learning methods. 

2.2 Identification Technique of Merging Systems. 

Galaxy systems were identified in the mock catalogue reproducing 
the FoF process as applied on real redshift catalogues ( [Merchan &| 
|Zandivar^P002[ (, |Merchan & Zandivarez| ( |2005[ ). 

Using the properties of the galaxies in these FoF systems we 
compute several features relevant to the problem, namely: 

(i) The DS test developed by jPressler & Shectman|(T988^ , uses 
the deviation of the local radial velocity, defined as the mean radial 
velocity of the closest n galaxies to each galaxy, from the global 
radial velocity in order to find substructures in clusters of galaxies. 
A global cluster value of the DS test is then obtained by summing 
up individual galaxy values. Following jPinkney et al.| ( |T996| l we se¬ 
lect those systems of galaxies with an occupancy Ngal > 30 galaxy 


With this set of features we test different machine learning 
algorithms such as the Logistic Regression ([Davison & Hinkley 


1997| Canty & Ripley|20T5), Su pport Vector Machines i Cortes & 
Vapni c|19^ Meyeret^J2014 1 and Random Forest (hereafter RF, 


Breiman[2001[ Liaw & Wiener 2002[l, provided by the R statistical 


programming language, with the aim of finding merging systems in 
the complete sample of galaxy clusters in our simulated catalogues. 

In order to measure its performance we run a standard cross 
validation test in 8 folds i.e. we divide the total sample in 8 in¬ 
dividual and independent subsets and train each machine learning 
algorithm with 7 of them in order to predict the dynamical status 
of the clusters of the remaining fold. As we know both, the under¬ 
lying and predicted dynamical status, we are able to compute the 
true positive rate (TPR), defined as the ratio between the number 
of merging clusters found in the final sample and the total num¬ 
ber of merging clusters that were in the test fold. In the same way 
we are able to compute the false positive rate (FPR), defined as the 
number of relaxed clusters classified as merging clusters divided by 
the number of relaxed clusters that the studied fold has. This infor¬ 
mation allow us to construct the Receiver Operating Characteristic 
(ROC) curve showed in Fig.[T](a) in which it can be seen that the 
best performance is obtained by the RF algorithm. 

For each cluster the RF computes a statistic that is related to 
the probability that the cluster is undergoing a merger, by building 
many decision trees from bootstrap training data, where the final 
classification is based on the average assignation of the ensemble 
of decision trees. Each tree is grown using randomly selected fea¬ 
tures from the training dataset previously described. We impose a 
threshold to the RF statistic of each cluster in order to classify the 
sample of merging clusters. In order to select an optimal thresh¬ 
old, we study in Fig. [T] (b) the impact of different values on sev¬ 
eral statistics of the sample of merging systems, namely the TPR, 
the FPR, the effectiveness (number of identified true-mergers di¬ 
vided by the total number of identified mergers) and the normal¬ 
ized length (number of identified mergers divided by the maximum 
number of identified mergers of the different thresholds). As can 
be seen, the classification threshold impacts on the performance of 
the RF classifier. Consequently, we select a threshold value of 0.3. 
This selection criterion guarantees a low false positive detection 
(high effectiveness) in the selected merging systems, but it should 
be recalled that we are only detecting just a fraction of the overall 
merging systems in the simulated catalogue as can be seen in Fig.[2 
(c). We note panel (d) of Fig.[T]will be analyze later in section 2.3. 


The RF implementation also allows us to assess the relative 
importance of the features as described in [Ehrlingerj ^2015[ l. In 
our case the most important features are the number of galaxy 
members, the p-value of the Shapiro-Wilk test and the Dressler- 
Shectman test. Nevertheless, it is worth noting that all the previ¬ 
ously described features are used in our RE implementation. 
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Figure 1. (a) ROC curve for logistic regression (red crosses), support vector machines (blue diamonds) and random forest (black dots) classifiers evaluated 
using 8 folds cross validation test on the SDSS mock training catalogues, (b) Effectiveness (black dots), true positive rate (red crosses), false positive rate 
(blue asterisks) and normalized length (green diamonds) as a function of the threshold imposed to the RF statistic (see text). Vertical line indicates the selected 
threshold value, (c) Effectiveness (black dots) and true positive rate (red crosses) as a function of the number of galaxy systems for the RF model, with the 
values at the selected threshold indicated by continuous lines, (d) Normalized projected distances Tnorm and velocity difference v„omi of the merging clusfer 
haloes recovered by our methodology (black points). False positive merging clusters are indicated by magenta squares. Merging clusters with missidentified 
substructures are indicated by red crosses. Vertical line indicates the selection cut introduced in order to avoid LOS contamination. 


2.3 Measured Properties of the Merging Haloes. 

The RF algorithm gives us a list of galaxies with high probability 
of belonging to the merging system. On the other hand, this method 
does not provide information about galaxy membership to the in¬ 
dividual substructures. Using a Mixture of Gaussian algorithms (R 
package mc/nvf, [Fraley et 5^|2012| l it is possible to cluster these 
member galaxies into the two merging substructures. Assigning 
galaxies we are able to compute the centre position (angular and 
redshift) of the substructures, its velocity dispersion (Gapper esti¬ 
mator), virial radius and, therefore, measure a dynamical mass. In 
order to estimate the associated errors we implemented a bootstrap 
technique. As shown in Fig.|^(a, b) the individual masses and the 
mass ratios are well recovered. 

The recovered geometries of the merging systems can be 
seen in Fig. [T] (d) in terms of the separation of the components 
in the line of sight (LOS) and in the plane of the sky. Specifi¬ 
cally in the x-axis we compute the normalized projected distance 
Gorm = ^i. 2 /(rviri + Girz). where di 2 is the angular separation be¬ 
tween both components of the merging system and Tviri and 
are the corresponding virial radii. In the y-axis we show the ve¬ 
locity difference Vnorm = \vi, 2 \l(o'i + 0 - 2 ) where Vi ,2 is the velocity 
distance between components and (Tj and cr 2 are the corresponding 
velocity dispersions of the substructures. Using this identification 
method we found 3 different cases. Relaxed clusters that we clas¬ 
sify as merging clusters (indicated by magenta squares), merging 
clusters in which we are unable to recover the real substructures 
(indicated by red crosses) and merging clusters in which we do 
recover the true substructures that are undergoing a merger (indi¬ 


cated by black dots). As can be seen the false positive cases (FPR 
~ 15%), indicated by magenta squares, are evenly distributed as a 
function of the parameter. The merging systems where we are 
not able to recover the real substructures are concentrated below a 
value of r„orm = 0.22. Such a selection cut will be used in order 
to report the properties of the substructures. Although we are able 
to find spatially coincident merging systems (r„orm < 0.22, that is 
systems merging along the LOS or systems with a small projected 
angular separation) we are not able to recover well the intervening 
substructures. 

Using this methodology we are able to define samples of 
merging systems with high levels of purity, low contamination and 
almost exact computation of the centre position of each compo¬ 
nent, as can be seen in Fig.[^(c,d). In panel (c) we show the angu¬ 
lar separation 0 between the actual and the recovered positions of 
each component, normalized to the actual virial radius. Similarly in 
panel (d) we show the differences Av between actual and measured 
radial velocities, divided by the real velocity dispersion. It should 
be recalled that our machine learning method is able to recover the 
correct substructures which belong to the major merger as identi¬ 
fied using the merger tree. 


3 APPLICATION TO LOW Z CLUSTERS. 

After testing our algorithm on simulated data, we applied our iden¬ 
tification algorithms to galaxy systems with more than 30 galaxy 
members identified in the update d ([Merchan & Zandivarez|2002| l 
catalogue, based on SDSS-DR7 jAbazajian et al.||2009^ data and 
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Figure 2. Recovered properties of merging galaxy systems identified in the simulated mock catalogue, (1) and (2) refers to the main and the merging subhaloes 
respectively: (a) Estimated mass bias (Mobs/Mrcai - 1) for the main and merging substructure. Color scale represent different mass ratios, (b) Probability 
distribution function for the real and observed mass ratio (M 2 /M 1 ). Panel (c) displays the cluster center separation on the sky plane (real and recovered) 
normalized to the real virial radius. Panel (d) shows the absolute value of the velocity separation (real and recovered) normalized to the real velocity dispersion. 
See text for further explanations. 


two samples of galaxy cluster measurements: WINGS ( |Cava et ^ 
|2009| > and HeCS ( |Rines et al.|2013| >. For the WINGS clusters we 
compute the g-r colour of the individual galaxies based on the ob¬ 
served b-v colour applying the formulas presented by the 2df Col¬ 
laboration based on the results of |Fukugita et al.| ( [l995| l. 

We report the following clusters as spatially coincident merg¬ 
ing systems candidates: A2593, A2199y, A2048y, A3266y, 
A3497 (full list will be available in the published version).Notice 
that we do not provide the substructure properties due to the LOS 
projection effect (as discussed in section \2.3\ , therefore minor 
mergers could be included in this sample ^Ma et ar]|2012^ . We 
display in Table [T] the properties of the two interacting structures 
of the merging clusters identified using our machine learning (RF) 
classification techniques. The errors for each property are the stan¬ 
dard deviation computed on a hundred of measurements obtained 
from the RF and the clustering realizations. Many of them are well 
known merging systems (indicated by /in the table), however it is 
important to emphasize that our method was able to find several 
new candidates and also to measure their properties. It should be 
noted that the decomposition of these structures is only indicative, 
because galaxy velocities are strongly affected by the gravitational 
attraction of the two halos. Therefore a tomographic reconstruction 
is necessary (including lensing and X-ray/SZ data) to recover the 
substructures accurately. 


3.1 The case of multiple major mergers. 

It is well known that there are some clusters that are the result of the 
merging of more than two systems, although they represent only 
a small fraction of the total sample (27 of 132 in 8 SDSS mock 
catalogues). In order to recover all the merging substructure we 


performed a mixture of more than two Gaussians. We found that 
our algorithm is only able to recover, with reliable properties, the 
two more important substructures. Hence the remaining structures 
may appear as contamination or not appear at all. We discuss an 
individual analysis of Abell 1758 as an example of the multiple 
major merger case. This cluster is known as a merger of four sub¬ 
structures, two in the north and two in the south ^Ragozzine et al .| 
|2012a^ . At a first iteration, our algorithm, classified correctly Abell 
1758 as a merging cluster, but failed to properly reconstruct the 
merging substructures. Considering this extra information we sep¬ 
arate the cluster in north and south components and perform a new 
analysis to each separately. We find that our algorithm was able 
to classify both as merging clusters and well recovers the merging 
substructure properties of both components (north and south). 


4 CONCLUSIONS AND FUTURE WORKS 

In this work we introduce a method aimed to detect merging sys¬ 
tems of galaxies in redshift surveys. We select a random forest 
algorithm between other machine learning algorithms, and use as 
features quantities derived from the galaxy redshift space informa¬ 
tion and from photometry (eg. colours). Our detection method was 
trained and calibrated using a sample of merging systems extracted 
from mock catalogues. By studying the merger trees we check that 
we do find the two substructures that experienced a major merger 
and recover their fundamental properties (positions and masses). 
We apply our techniques to a sample of systems of galaxies identi¬ 
fied in SDSS-DR7, WINGS and HeCS. 

The resulting merging system sample in which we are able to 
recover the merging substructures comprises 12, 4 and 16 systems 
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respectively. Additionally, we report 29 spatially coincident merg¬ 
ing system candidates. Several of this systems where previously 
reported by other authors as interacting systems of galaxies. We 
also report for the first time 40 new nearby candidates as merging 
systems, that were previously overlooked. 

We emphasize that our method detects in a reliable way the 
merging systems candidates and substructures properties, but also 
wish to note that, in the case of multiple mergers, some merging 
substructures may be joined by our algorithm and hence are re¬ 
ported as one component causing a possible biase in some mea¬ 
sured properties. 

The kinematical reconstructions (see Fig. [^and Table [g will 
be corroborated using tomographic techniques including X-ray and 
weak lensing information jGonzalez et al.|20T5l l in forthcoming pa¬ 
pers. We also plan in further works to train our algorithm with light 
cone mocks in order to apply our technique to high redshift cata¬ 
logues (CLASH-VLT (B iviano et al. pOT^, FRONTIERS (Ebelmg] 
|et al.|2014T l, EDiCs jMilvang-Jensen et al. 2008| l, DESI, etc.) and to 
study any environmental dependence of galaxy properties (eg. star 
formation rate, stellar mass, morphology) at different stages in the 
merging process. Diverse studies could be performed with a sample 
of merging systems like the ones presented in this paper. Using the 
Bayesian reconstruction techniques presented by |Dawson| ( |2013^ it 
is possible to recover the 3D information of the merger. In a forth¬ 
coming paper we will use such information in order to explore the 
implications for the properties of the DM particle. A web inter¬ 
face implementing these methods (The MeSsI Algorithm) is freely 
available at http://200.16.29.98/martin/merclust 
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Table 1. Here we present the low redshift merging cluster sample, including the properties of the two main substructures identified by our algorithm. In column 
1 we present the name of the cluster, from column 2 to 5 we present the estimated mass and the position of the main substructures and from column 6 to 9 we 
present the estimated mass and the position of the other substructure, finally in the last column we list previous work on each cluster. Clusters that have been 


previously reported as merging systems are indicated with /. References: 2 

Wen & HanH2013|, 3 

Einasto et al.H2012), 4|Cohen et al.H2014), 6|Abdullah 


201 1^, 7|Rines & DiaferioH2006(, 8|Rhee et al. 

jl991k 9|Parekh et al.| 

2015), 16|Ramella et a 

.1(2007), 21|Wang et al.|(2010), 22|Johnston-Hollitt et al. 

j2008 

, 23|Boschin et al.H2012), 26|Ragozzine et al.| 

2012b), 27|Durret et aLH2011), 29|Smith & TaylorH2008), 30|Komgut et al. j^Olll. Full version of the 


Name 

Ml [lO’^'Mo] 

RAiP] 

DECi [°] 

Zl 

M2 [IO'^'Mo] 

RA2[“] 

DECiH 

Z2 

References 

Abell 1424 

4.9 

179.38 

5.08 

0.0760 

5.1 

179.19 

5.01 

0.0746 

3,6.7,8 


±2.3 

±0.09 

±0.02 

±0.0004 

±1.4 

±0.1 

±0.04 

±0.0005 


Abell 2142 

18.3 

239.61 

27.23 

0.0901 

11.3 

239.33 

27.5 

0.0893 

2 

y 

±0.6 

±0.005 

±0.005 

±0.0004 

±1.8 

±0.005 

±0.005 

±0.0001 


Abell 3158 

37.24 

55.75 

-53.63 

0.0633 

4.6 

55.37 

-53.48 

0.0622 

4,9,16,21,22 

y 

±1.5 

±0.07 

±0.004 

±0.0001 

±0.2 

±0.007 

±0.001 

±0.0001 


Abell 2382 

77.7 

327.90 

-15.66 

0.0676 

6.12 

328.167 

-15.62 

0.0642 



±10.2 

±0.006 

±0.006 

±0.0003 

±1.1 

±0.003 

±0.01 

±0.0002 


Abell 1758N 
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203.07 

50.59 

0.2768 

29.1 

203.25 

50.57 

0.2783 

2,23,26,27 
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±9 

±0.02 

±0.01 

±0.0002 

±15.8 

±0.02 

±0.03 

±0.0007 


Abell 1835 

105 

210.25 

2.87 

0.2516 

17 

210.29 

2.75 

0.2479 

2,29,30 


±14 

±0.01 

±0.01 

±0.0006 

±25 

±0.09 

±0.04 

±0.002 








Figure 3. Galaxy Angular distribution of some merging systems. From left to right we show on columns: two simulated systems, two clusters from the SDSS, 
two WINGs clusters and two HECs clusters. Member galaxies of both substructures are showed in red and blue dots while in black lines are plotted the 
iso-density contours. The ellipses indicate the 1, 2 and 3 sigma error of the identified substructure (on blue or red colours depending on the component), with 
solid, dashed and dotted line, respectively. For comparison, in green dots are plotted the angular positions of the X-ray sources, and in magenta dots are plotted 
the angular positions of the substructures identified by other authors. For the simulated merging systems we also show in dashed line the iso-density of the 
true substructures that are colliding and in red and blue lines de iso-density contours for the identified substructures. 
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